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ABSTRACT 


An  approximate  analytic  solution  has  been  developed  for  the 
ablated  volume  and  shape  of  a  metal  sphere  at  hypersonic  speeds 
in  the  early  stages  of  aerodynamic  ablation.  The  process  of  abla¬ 
tion,  due  to  aerodynamic  heating,  is  a  very  complicated  phenom¬ 
enon  because  of  the  possible  existence  of  two  layers  -  gas  and 
liquid.  If  the  metal  surface  is  molten,  then  a  liquid  layer  of  mol¬ 
ten  metal  is  formed  between  the  boundary  layer  of  gas  and  the 
solid  surface.  The  interdependent  factors  flight  time,  flight  con¬ 
ditions,  and  body  shape  cause  a  variation  of  aerodynamic  charac¬ 
teristics  as  a  function  of  the  trajectory.  Furthermore,  the  aero¬ 
dynamic  ablation  is  responsible  for  a  change  in  the  shape  of  the 
body.  Two  cases  of  aerodynamic  ablation  were  investigated.  The 
first  case  considered  the  melting  of  body  material  alone,  and 
the  second  considered  melting  with  partial  evaporation. 

The  analysis  given  in  this  report  is  useful  to  determine  the 
flight  conditions  under  i^ch  the  metal  sphere  may  travel  with 
little  or  no  ablation.  With  the  assumption  that  the  flow  over  the 
sphere  is  laminar,  Newton’s  impact  theory  is  applicable  and  the 
gas  does  not  dissociate;  the  total  ablated  volume  and  shape  of 
the  sphere  can  be  computed  as  a  function  of  flight  trajectory  by 
using  the  analysis  developed  in  this  report. 
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ANALYSIS  OF  THE  AERODYNAMIC 
ABLATION  OF  A  METAL  SPHERE 


INTRODUCTION 

The  natural  phenomenon  of  ablation  can  be  observed  at  night  when  meteoroids  burn 
because  of  aerocfynamic  heating  and  leave  luminous  trails.  The  problem  of  aerotfynamic 
ablation  is  an  important  scientific  and  technical  subject  of  many  aerospace  programs. 

During  hypersonic  flight  in  an  atmospheric  environment,  excessive  heat  is  generated. 
Structures  intended  for  high-speed  flight  must  be  protected  from  this  intense  heating. 
Ablation  is  a  beat-  and  mass-transfer  phenomenon,  involving  a  physical  and  chemical 
process  which  protects  the  structirre  from  disintegration  by  effective  removal  of  heat 
from  the  body  surface.  When  the  wall  temperature  of  a  metallic  body  reaches  the  melting 
point,  the  surface  material  begins  to  melt,  and  a  liquid  layer  is  formed.  The  melt  is 
swept  downstream  by  aerodynamic  pressure  and  shear  forces.  If  the  temperature  of  the 
surface  increases  further,  the  surface  material  evaporates  and  the  metallic  gas  pene¬ 
trates  into  the  surrounding  gas  layer,  reducing  the  heat  flux  into  the  solid  body. 

A  sxirvey  of  available  references  on  aerodynamic  ablation  shows  that  the  problem  of 
ablation  has  be^  treated  for  the  most  part  analytically,  since  experimental  results  are 
scarce.  Goodman  (1)  developed  an  analysis  of  the  liquid  layer  which  takes  into  account 
the  heat  penetration  into  the  solid  body.  Bethe  and  Adams  (2)  and  Hidalgo  (3)  analyzed 
the  ablation  of  glassy  materials,  and  Scala  (4),  Roberts  (5),  and  Goodman  (6)  solved  the 
stagnation  point  ablation  problem.  Goodman  (7)  found  a  particular  solution  for  the  abla¬ 
tion  of  a  flat  plate  with  a  constant  heat  flux,  and  Roberts  (8)  formulated  an  engineering 
solution  for  the  sublimation  of  a  blunt  body. 

A  sphere  is  widely  used  for  aerodynamic  experiments,  and  its  blimted  nose  absorbs 
a  high  heat  flux.  Some  experimental  data  on  ablation  of  metal  spheres  are  available  at 
NHL.  A  theoretical  analysis  of  the  ablation  of  a  metal  sphere  is  desired.  Depending 
iq)on  flight  conditions,  three  cases  of  ablation  may  occur:  ablation  due  to  melting  only, 
ablation  due  to  melting  and  partial  evaporation,  and  ablation  due  to  sublimation  only. 
Roberts  (8)  solved  the  sublimation-only  case  with  an  engineering  solution. 

After  examination  of  the  available  analytical  and  experimental  references  on  the 
problem,  an  attempt  is  made  in  this  report  to  predict  the  total  volume  ablated  and  the 
new  shape  of  the  sphere,  due  to  melting  only,  and  due  to  melting  and  partial  evaporation. 
The  solution  presented  is  limited  to  the  early  stages  of  ablation  of  a  metal  sphere.  For 
the  analysis,  it  is  assumed  that  the  gas  and  liquid  flow  are  laminar,  and  that  the  gas  flow 
is  Newtonian,  Although,  because  of  the  complicated  process  of  ablation  which  depends  on 
the  cliaracteristics  of  flow,  time,  and  the  shape  of  the  body,  only  the  initial  condition  of 
ablation  is  analyzed,  the  developed  analytical  method  is  useful  to  designers  and  research 
scientists  because  the  flight  conditions  for  zero  ablation  can  be  determined  by  this 
method.  For  example,  the  ambient  pressure  for  flight  with  zero  ablation  or  the  time  at 
which  ablation  starts  can  be  computed.  For  such  calculations  the  velocity  and  trajectory 
are  assumed  to  be  known  as  a  function  of  time.  Therefore,  a  steady-flow  analysis  can  be 


Note:  Dr.  Chang  is  a  faculty  member  of  the  Department  of  Mechanical  Engineering,  The 
Catholic  University  of  America.  He  is  a  part  time  employee  of  the  Naval  Research 
Laboratory. 
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carried  out,  because  for  small  time  intervals,  the  characteristics  of  flow  and  the  shape 
of  the  body  may  be  assumed  independent  of  time.  It  is  further  assumed  that  the  gas  is 
nondissociative.  The  reasons  for  this  assumption  are  as  follorws:  Dissociation  occurs  if 
the  temperature  is  high  and  high-energy  molecules  collide.  For  molecules  of  oxygen  and 
nitrogen,  of  which  air  consists,  the  dissociation  occurs  for  oxygen  at  3000‘’K  and  nitrogen 
at  5000'’K,  when  these  molecules  collide  with  others  of  the  lowest  minimum  relative  ki¬ 
netic  energy  D  (9).  The  values  of  DAT  for  these  elements  are  59,300/t  for  oxygen  and 
113,100/t  for  nitrogen,  where  k  Is  Boltzmann's  constant  (l.38xl0'*6  erg/molecule-®K). 
However,  molecules  with  such  high  relative  kinetic  energies  are  rare,  therefore,  in  order 
for  dissociation  to  occur,  a  considerable  time  is  required.  This  is  especially  true  for 
nitrogen  (9),  This  report  deals  with  short-time  ablation  in  the  order  of  milliseconds; 
therefore,  dissociation  of  the  gas  has  been  neglected.  This  anatysis  is  applicable  for 
very  thin  liquid  layers  in  comparison  to  the  radius  of  the  sphere  and  in  the  range  of  Mach 
numbers  M  >  6,  where  M  refers  to  the  free  stream  in  front  of  the  shock  wave. 


THE  EQUATION  OF  THE  LIQUID  LAYER 


h. 


From  Fig.  1,  it  is  seen  that  the  volume  of  ablated  material  depends  on  the  distance 
However,  because  h  =  s  -  5,  in  order  to  determine  h,  it  is  sufficient  to  compute  s 

and  ^ .  Therefore,  for  the  determination  of 
the  thickness  of  Hie  liquid  layer  the 
equation  of  the  liquid  layer  has  to  be 
derived. 


Fig.  1  -  Sphere  in  hypersonic  flight 


When  melting  occurs,  a  metallic  body 
has  a  definite  sur&ce  of  demarcation  be¬ 
tween  the  metal  liquid  and  the  solid,  in 
contrast  to  a  glassy  heat- shield  material 
which  possesses  no  distinct  solid  state. 
The  analytic  solution  for  the  total  volume 
of  material  ablated  and  the  final  shape  of 
the  sphere  wiU  be  presented  for  steady- 
state  ablation  and  laminar  flow. 

The  equations  governing  the  liquid 
layer  around  a  ^here  are  essentially  the 
same  as  those  for  incompressible  fluid 
flow  around  a  body  of  revolution.  Initially 
ablation  takes  place  as  Fig.  1  indicates, 
where  the  applicable  coordinate  system  is 
also  shown.  The  coordinate  x  along  the 
surface  of  the  sphere  is  measured  from 
the  stagnation  point,  and  the  coordinate  -y 
is  perpendicular  to  x.  At  the  early  stage 
of  ablation,  foe  loss  of  mass  is  small  and 
foe  small  arc  length  x  may  be  replaced  by 
a  straight  segment  on  the  gas/Iiquld  inter¬ 
face.  For  foe  computation  the  coordinate, 
y  is  assumed  to  be  parallel  to  h  and  s. 

For  steady-stdte  incompressible  flow, 
the  governing  equations  are 


Continuity; 


rv) 

dy 
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3 


Energy;  pc 


.  ■  3u  3u  3  /  3u\ 

Momentum:  pu  ^  ^ 


where  I  Is  the  radiation  heat  flux.*  Rdberts  (10)  found  that  for  a  ballistic  vehicle  the 
heat  lost  by  radiation  is  negligible  compared  with  the  heat  disposed  of  by  ablation.  For 
simplicity,  It  Is  assumed  that  the  dynamic  viscosity  p  is  constant.  Because  the  viscos¬ 
ity  of  the  liquid  layer  is  high  and  the  inertia  forces  of  the  same  layer  are  negligible,  the 
terms  of  puf3u/3x)  and  pvf3u/3y)  in  the  momentum  equation  as  well  as  u(3p/3x)  and 
/i(3u/3y)2  in  the  energy  equation  are  neglected.  The  gradient  of  temperature  3t/3x  is 
small  compared  to  3T/3y,  and  since  u  and  v  are  also  small,  only  the  term  3T/3y  is 
considered  in  the  energy  equation.  Then  the  governing  equations  become 


n) 


Equation  (2)  shows  that  the  heat  flux  through  the  liquid  layer  is  a  constant;  therefore 


3*u  3p 
^2  3x 


From  Eq.  (1)  for  the  liquid  layer 


Because  the  liquid  layer  is  constantly  swept  downstream  and  is  kept  thin,  r,  the  radial 
coordinate  from  the  center  line  to  the  surface,  is  regarded  as  constant.  Hence 

I  dy  =  -r[v(s)  -  v(h)] 


-  r 

3x  J 


(ru)  dy  =  f  dy  +  u(s)r  ^  -  ru(h)  ^ 

\  «  3x 


Since  the  velocity  on  the  surface  of  the  body  is  zero,  it  follows  that  uf  s)  =  0 ,  and  since 
in  the  outer  layer  the  vertical  velocity  v  consists  of  local  and  convective  velocity  com¬ 
ponents. 


*A  list  of  definitions  of  the  symbols  follows  the  main  body  of  this  report. 
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^  I’fru)  dy 


f’ 

-  ;  a* 


3h 

<ly  +  rv{h)  -  r  -5^ 


From  Eq.  (4) 


A  I  (ru)dy+r[v(.)-|t]=  0. 


If  there  exists  a  convective  velocity  caused  by  the  change  of  density  between  liquid  and 
solid,  then  the  normal  velocity  component  in  the  liquid  side  of  the  liquid/solid  interface 
is  given  (11)  by 


/  ^»ol \  Bs 


Denoting  5  =  s  -  h,  Eq.  (5)  is  written  as 


,  ^  ^  3^  '°.ol  3s 

f  ru)  dy  +  r  —  =  r  -  —  . 


Distribution  of  the  Tangential  Velocity 

In  order  to  solve  Eq.  (7),  it  is  necessary  to  determine  the  tangential  velocity  u.  The 
quantity  u  is  obtained  by  using  the  momentum  equation 


3x  ■  3y2  • 


By  integrating  and  using  the  boundary  condition 


Bu 

^3y  “ 


at  y  =  h ,  it  follows  that 


By  integrating  once  more  with  respect  to  y  and  using  the  boundary  condition  u(s)  =  0  at 
y  =  s,  the  velocity  component  in  the  x-direcUon  becomes 


Governing  Equation  of  the  Liquid  Layer 

The  thickness  of  the  liquid  layer,  can  be  determined  by  using  the  governing  equa¬ 
tion  of  the  liquid  layer. 
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Substituting  Eq.  (8)  into  Eq.  (7), 


3x 


♦  r 


Ik 

Bt 


^3 

ST  ' 


(9) 


IMs  is  the  same  equation  as  derived  by  Goodman  (6).  The  right-hand  side  of  Eq,  (9)  ex¬ 
presses  the  rate  of  melt  production,  which  is  a  function  of  external  heat  flux  and  the  rate 
c€  heat  penetration  into  the  solid.  The  effect  of  heat  penetration  into  the  solid  due  to 
aerodynamic  ablation  is  given  by  Goodman  (1)  as 


Bs  r  2  1 

^*oi  ”  at  ■  ■ 

Then  Eq.  (9)  is  reduced  to 


1  B  r 

0* 

k 

1 

IsJi  I 

rt)  I 

t 

A 

1 — 

1--] 

3x 

[  3  Bx  ^  +  2  ^ 

,J  ■  ^ 

e) 

wbere  0  is  given  implicitly  in  terms  of  F  and  v  by 


(10) 


(11) 


F 


-  2  +  2(  1  +  »-)  In 


2(l  +  v)  -  0 
2v 


F 


■t/(p  ,Hk  ,T.  ) 
sol  sol  is' 


(12) 


and 


Patting 


1  -  |  =  Gjfx.t) 


Bq.  (11)  becomes  (for  steady  conditions,  where  B^/Bt  =  O) 


1  !p  Ir 

3  3x  ^  2  ^ 


r 


Equation  (14)  may  be  regarded  as  the  final  equation  of  the  liquid  layer. 


(13) 


(14) 


AERODYNAMIC  ABLATION  OF  A  METAL  SPHERE 

In  the  case  of  ablation  due  to  melting  only,  the  liquid  viscosity  is  much  larger  than 
tttat  of  the  gas;  therefore  the  tangential  veloci^  of  the  liquid  is  small  compared  with  the 
gas  velocity.  Hence  the  existence  of  a  liquid  layer  may  be  ignored  for  the  evaluation  of 
property  values  of  the  gas  layer,  at  the  Initial  condition  of  ablation.  The  pressure  dis¬ 
tribution  of  gas  layer  over  the  liquid  layer  is  assxuned  to  be  equal  to  that  over  the  solid 
surface.  However,  if  evaporation  occurs  in  addition  to  melting,  then  the  evaporated 
ntetal  gas  penetrates  into  the  gas  layer,  and  the  gas  layer  characteristics  will  change. 
Ibis  effect  is  given  by  shear  stress  ratio 


swJ,  using  Reynolds  analogy ^ 
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where  the  subscript  i  refers  to  the  gas/llquid  Interface  and  the  subscript  o  to  the  non¬ 
ablating  body.  The  effect  of  gas-molecule  dissociation  due  to  aerodynamic  heating  at 
high  speeds  Is  neglected  In  the  analysis.  This  assumption  may  hold  for  the  short  atmos¬ 
pheric  flight  path  of  20  meters  and  the  maximum  Mach  number  of  about  18  of  the  NRL 
test  facility,  because  nitrogen  dissociation  is  not  expected  to  be  serious  if  it  occurs  (12), 
The  gas  and  liquid  layers  are  assumed  to  be  laminar. 


Ablation  Due  to  Melting  Only 

In  order  to  solve  Eq.  (14)  of  the  liquid  layer,  the  liquid-flow  properties  contained  in 
this  equation  are  determined  through  the  gas-flow  properties  around  the  sphere.  Since 
the  thin  liquid  layer  is  essentially  a  boundary  layer,  and  since  the  static  pressure  of  the 
liquid  is  equal  to  that  of  the  ambient  gas,  the  pressure  distribution  of  the  liquid  layer  can 
be  formed  by  that  of  the  gas  around  the  sphere.  The  Newtonian  impact  theory  at  hyper¬ 
sonic  velocities  yields  the  pressure  coefficient  around  the  sphere  as  =  2  sin^  S, 
where  8  is  the  angle  between  the  tangent  to  its  local  surface  of  the  body  and  the  stream 
direction.  Lees  (13)  proposed  a  modified  Newtonian  theory,  which  consists  in  scaling 
down  this  result  so  as  to  be  exact  at  the  stagnation  point  of  a  blunt  bo^,  at  which  point 
the  correct  value  is  known.  Then  the  modified  Newtonian  relation  is. 


This  simple  equation  is  in  a  very  good  agreement  with  experimental  data  in  the  region 
0  <  x/S  <  1.4.  Hence,  the  pressure  distribution  around  the  sphere  is  determined  from 
Eq.  (16): 

P  -  P„  =  (P,~P„)  cos*  (x/R)  . 

Then 


3p 

ax 


fPs*-P„> 

R 


Shear  Stress  in  the  Liquid  Layer 

It  is  necessary  to  evaluate  the  shear  stress  in  the  liquid  layer  in  Eq.  (14) .  This 
value  can  be,  however,  expressed  by  the  shear  stress  of  the  gas  layer.  The  shear  stress 
in  the  gas  boundary  layer  can  be  computed  if  the  thin  liquid  layer  is  ignored  and  the  ve¬ 
locity  profile  is  assumed  to  be  the  shape  of  the  dashed  curve  in  Fig.  2.  However,  such  a 
velocity  distribution  deviates  from  the  real  velocity  distribution  in  the  gas  boundary 
layer,  so  the  real  shear  stress  is  given  by 


where  and  tj,  are  real  and  ideal  shear  stresses  in  the  gas  and  correspond  to  the 
shear  stresses  of  the  velocity  distributions  of  the  gas  boundary  layer  shown  by  full  and 
broken  curves  respectively  in  Fig,  2.  K  is  a  factor  to  be  determined.  Van  Driest  (14) 
shows  that  for  laminar  gas  flow,  for  a  sphere. 
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Accounting  for  heat  transfer,  may  be  writ¬ 
ten  as 

’’u  =  (l.  32/^11,  X 

where  the  superscript  *  refers  to  the  property 
values  evaluated  at  the  reference  temperature 
(15) 

T’/T,  =  1+0.035  m/  +  0.45  [(Tyt.)  -  l]  .  (1"^^) 


GAS  BOUMDARY  LAYER 


LIOOID  LAYER 


Fig.  2  -  Velocity  distribution 


The  Mach  number  M,  along  the  sphere  surface  behind  the  detached  shock  is  approximated 
(16)  by 


"  2  R 


Using  the  isentropic  flow  relationship  in  the  region  between  the  shock  wave  and  the  body 
surface. 


Then  for  a  diatomic  gas,  y  =  1.4  and 


3  27  3  x3  A  27  x2\ 

“ej,  8  r3  ■  40  gaj 


and  for  diatomic  gases  at  temperatures  higher  than  2000“K, 


Therefore,  from  Eq.  (17) 


For  laminar  flow,  the  ratio  of  shear  stresses  may  be  written  as 


^liq  _  ^liq 

T  ~  U 

^gas 


where 
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The  dynarr;tc  viscosity  of  molten  metal  is  given  by 


where  A  and  b  are  constants.  For  aluminum  A  =  705°K  and  b  =  0.474,  and  for  Iron 
A  =  1690^K  and  B  =  0.854  (17), 

Since 


For  a  diatomic  gas 


Differential  Equation  of  the  Liquid  Layer  and  Determination 
of  the  Liquid  Layer  Thickness 

If  in  addition  to  the  information  of  the  preceding  section  the  heat  flux  over  the  sphere 
is  known,  then  the  differential  equation  of  the  liquid  layer  can  be  derived.  The  heat  flux 
around  the  hemisphere  is  given  (18)  by* 

<1  =  1.  (I)]  I  l}/  [f.  (f )  *  (f )]‘^ 

=  -.={!) 

where  ^  =  l-fp'/l?,)  and  ?*  is  the  static  pressure  at  the  shoulder,  and  where 
♦The  power  1/2  In  Eq.  (22)  was  erroneovisdy  omitted  in  Ref.  18. 
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(f)  =  -  -  ii)  (I)  [i  *  (t)]  •  f  ■‘"’(k)*  5  •‘■'*(1)}  • 


For  M  >  6,  /3  =  0,96, 


The  stagnation  heat  flux  is  computed  from  Ref,  19: 

0,763  (a^y 

■  p,o.6\u„y  \p„o.Dj 

Taking  =  1,1  (u„/D)  for  a  sphere. 


/°a>“<DCp(T,  -T,). 


■  Pr0.6(p„-„D)  (pj  r 

After  substitution  of  Eqs.  (16),  (21),  and  (22),  Eq,  (14)  becomes 

1  fd?;  ~  1 


where 


^2  =  (P.-Po.) 


2cos(f)sin(j)] 

^2  =  (P.-Po.) 

^3  =  f  [cos  (f)  G/x)  +  R  sin  (|)^^’] 

^4  =  aR  sin  j  GjY  x)  . 


In  order  to  estimate  the  magnitude  of  each  term,  a  numerical  computation  is  carried 
out  at  the  sonic  point  (see  Appendix  A).  The  following  assumptions  are  made:  a  4-cm- 
diameter  aluminum  sphere  travels  at  =  15  for  0,003  sec  in  a  standard  atmosphere; 
the  thickness  of  the  liquid  layer  is  C  =  0,001  cm.  The  results  show  the  following  orders 
of  magnitudes: 


f,  =  0(10’)  .  =  O(105)  . 


Since  the  magnitudes  of  the  terms  and  are  about  10%  of  the  terms  and 
the  terms  ^3  and  are  neglected.  Physically  S3  and  ^4  correspond  to  shear  resist¬ 
ance;  on  the  other  hand  and  ^ j  express  the  pressure  distribution  around  the  sphere. 
From  these  numerical  examples  it  can  be  seen  that  the  streamwise  decreasing  pressure 
plays  a  much  greater  role  in  sweeping  away  the  molten  layer  than  does  the  shear  force. 
Equation  (24)  is  then  reduced  to 
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^  +  ^(x)  ^  =  iVfx,  t)  5  ^  (25) 

where 


and 


^(x)  =  ^,/^3 


^x,  t) 


R/*nqq. 

^UqH 


sin(x/R)  r(x/R) 
^2 


Gjfx.t)  . 


By  taking  time  k  as  constant,  i.e.,  by  fixing  t,  v  becomes  a  function  of  x  only,  and  Eq, 
(25)  can  be  integrated.  Its  solution  is  given  by 

The  constant  c,  may  be  determined  by  the  boundary  and  initial  conditions,  C  =  0  at 
X  =  0  and  t  =  0,  The  thickness  of  the  liquid  layer  can  be  calculated  by  Eq.  (26)  at  con¬ 
stant  time  t .  Tine  shape  of  an  ablated  sphere  is  determined  by 

h  =  s  -  5 

where  s  is  obtained  by  integrating  Eq.  (10)  as 


1/3 

(26) 


s 


/  rdt\ 


q,g(x/R) 


+  c, 


(27) 


The  boundary  and  initial  conditions  are  at  t  =  0,  for  all  x,  s  =  0.  These  conditions  can 
be  satisfied  if  c^x)  is  constant.  Therefore 


Then 


t»3(x)  +  Cj  =  C3  . 


s 


q.gCx/R) 

^liqH 


+  C, 


The  volume  of  ablated  material  V*  is  at  the  fixed  time  t 


▼*  =  2'n’  J*”  hfj  dx  =  (s-  C)  I^R  -  ^ (^}*^** 

Thus  far  C,  h,  and  v*  have  been  expressed  by  stagnation  conditions.  However,  the  stag¬ 
nation  properties  are  determined  by  the  free  stream  values  in  front  of  the  shock  and  by 
normal  shock  relations.  Hence,  the  total  ablated  volume  and  shape  of  the  sphere  are 
evaluated  by  tbe  trajectory  at  any  specified  time  t . 

The  physical  meaning  of  the  equations  derived  is  summarized  as  follows:  If  the  ve¬ 
locity  of  the  sphere  is  zero,  but  a  heat  flux  exists,  the  molten  layer  s  increases  with 
time  t  and  the  thickness  C  is  equal  to  s ,  since  no  ablation  occurs  under  such  conditions. 
However,  if  the  velocity  of  the  sphere  is  not  zero,  then  the  molten  layer  is  swept  away 
because  of  the  pressure  decrease  in  the  streamwise  direction,  and  h  increases  with  time. 
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In  other  words  C  will  decrease.  As  Eq.  (28)  shows,  with  decreasing  thickness  of  the 
liquid  layer,  the  ablated  volume  increases  accordingly.  The  thickness  ^  can  be  deter¬ 
mined  by  solving  Eq.  (24)  as  a  closed  form  of  Eq.  (26)  as  a  function  of  x,  if  the  time  is 
fixed.  Therefore,  it  is  possible  to  calculate  numerically,  using  Eq.  (28),  the  total  ablated 
volume  at  any  fixed  time  of  the  trajectory.  The  determination  of  the  ablated  mass  at  a 
fixed  time  is  in  practice  very  important,  for  example,  for  problems  of  heat  transfer,  de¬ 
sign  of  high  speed  vehicles.  Impact  of  bodies  traveling  at  hypersonic  speed,  and  intensi¬ 
ties  of  luminous  trails  and  ionization  behind  bodies. 


Ablation  Due  to  Melting  with  Partial  Evaporation 

The  occurrence  of  evaporation  may  be  determined  by  the  wall  pressure  or  conduc¬ 
tion  parameter.  For  an  ice  sphere  immersed  in  high-speed  air,  it  is  known  that  for  a 
given  recovery  temperature  there  exists  a  value  of  wall  pressure  at  which  the  wall  tem¬ 
perature  is  just  at  the  melting  point  of  the  ice.  For  wall  pressures  above  this  critical 
pressure,  the  ice  melts,  but  if  the  pressure  at  the  wall  is  less  than  the  critical  pressure, 
then  the  ice  does  not  melt  but  evaporates  (20).  The  dimensionless  conduction  parameter 
V  is  given  (8)  by 

[ra/{CoAS)]^^^ 

where  the  subscript  b  refers  to  body.  If  v  >  2x10^,  then  no  sublimation  occurs.  Bethe 
and  Adams  (2)  worked  out  the  theory  for  aerod|ynamic  ablation  with  partial  evaporation, 
and  their  formulations  are  applied  here. 

For  partial  evaporation,  the  ratio  of  shear  stress  or  heat  flux  at  the  gas/liquid 
interface  to  that  for  a  nonablating  body,  Eq.  (15),  is  an  important  parameter.  In  a  strong 
evaporation 


- 

h„  +  0.68Mi  (hi.-hi.) 


The  metal  vapor  pressure  at  the  gas/liquid  interface  is 


P  ,  =■ 

Vi 


-  0.68  M, 


The  quantity  Vj ,  the  normal  velocity  of  evaporated  metal  perpendicular  to  the  sphere 
surface  at  the  gas/liquid  interface,  is  determined  as  follows:  Knee,  in  addition  to  Eq. 
(29), 

su  1  -  0  68  M®'**  ^UqVi/nsgfx/R) 

the  two  equations  for  >/»  are  equated  to  give 


1  -  +  (hi  -hi  jj 

0.68  Mj'**  (’'l.'Nj) 
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Substituting  Eq.  (32)  into  Eq,  (30), 


= 


0,68  M, 


H  fl - 

.  0, 26 

+  0.68  Mj  (hi  -  hij 


-  0.68  M, 


Thus  p^j  Is  a  function  of  x. 


With  this  Information,  Eq.  (14)  is  reduced  to  the  following  form  if  small  terms  are 
neglected: 


d£_  __  sin  fx/R)  Misq'/'  q.gfx/R)  /2  \ 

d,  ♦  ^1  5  »  — j;;;® - (e-  ij 


where 


J  -  •  /^  \  ^vi 

and  a  bar  over  a  symbol  refers  to  ablation  with  partial  evaporation.  Then  it  follows  that 


where 


^  =  A  cot  (iU  i  ^ 

^2  3R  [kJ  3  3x2  /  Bx 


_  '/'qsB(x/R)  /J 

^  “  r»  O  in  ^  I  n 


^liqM  dp, 


For  a  fixed  time,  the  solution  of  Eq,  (34)  is 


VK)  /2  \ 

1“  («-■)• 


The  aerodynamic  heat  flux  is  balanced  by  the  energy  of  evaporation  and  tbe  heat  caoacitv 
of  the  liquid  and  solid  phases.  Thus 

*^«h*orptlon  ^  ^•vaporation  ~  ^liq  I'yl  '  (36) 
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Using  Eq.  (32), 


q  “  o  y 

•vflporation  *  '' 


q,  rC x/R)  hy  ^  hy 

0  ^  0  2/i 

).68  Mj'  ^hi  +  0.68  Mj'  ^hj  -hj^j 


The  absorbed  heat  can  be  computed  by 


'ttbsorpt 


=  d/Q  —  Q  ,  • 

ion  ^  ^or-iporation 


The  mass  loss  is  caused  by  evaporation  and  absorption,  so 


M  ^  S 

•vaporoti on  absorption 


ovaporat  ton 


abaorpt 


^evaporation  A  ^  f  dt\ 

tion  "  ^71  ii  7^  (*  ”  ^  J  ^  j 

^liq  evaporation  \  “^  / 

^absorption  /  ^  Tdt  \  . 

^~b  r  \  ^  J  ^  j 

iiq  absorption  \  / 


Hence,  the  shape  of  an  ablated  sphere  with  partial  evaporation  is  determined  by 


h  =  a  -  C 


^liq  /VH, 


evaporation  ^absorpt 

l""'"  ~r  U 

evaporation  absorpt 


ion  IT 

—  +  cJ  -  5. 

ioti  / 


The  total  ablated  volume  v*  is 


f*  =  27r  ^  h  Tj  dx  =  2w  (s  -  C)  j  (f")  dx  . 


Hence,  the  total  ablated  volume  and  shape  can  be  computed  based  upon  the  trajectory  at 
any  specified  time. 


DISCUSSION 

Aerodynamic  ablation  is  a  very  complex  phenomenon,  because  the  change  of  body 
shape  due  to  ablation  also  causes  changes  in  gas  and  liquid  flow  characteristics  along  the 
bo^.  The  approximate  analysis  presented  here  is  for  steady- state  aerodynamic  ablation 
of  a  sphere,  assuming  that  the  flow  characteristics  are  those  on  the  sphere  and  neglecting 
the  change  of  shape.  Hence,  this  analysis  is  applicable  for  a  flight  condition  with  small 
ablation.  This  method  is  especially  useful  to  determine  the  flight  conditions  imder  which 
the  ^here  can  travel  without  any  ablation.  However,  if  considerable  ablation  occurs,  the 
front  part  of  the  sphere  deforms  to  a  cone  shape  (Fig.  3).  No  attempt  has  been  made  to 
confirm  the  deformation  to  a  cone  by  analysis.  The  computation  of  change  of  shape  due  to 
ablation  will  be  possible  only  when  a  high-^>€^  computer  is  used.  If  the  moving  bounda¬ 
ries  for  melting  and  subliming  are  considered,  theh  the  problem  is  more  complex,  and 
Goodman’s  "heat-balance  integral"  (21)  may  be  useful  to  solve  the  problem  of  ablation 
with  changing  shape. 
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Fig,  3  -  Recovered  ablated 
aluminum  sphere  (original  di¬ 
ameter,  0.635  cm;  initial  ve¬ 
locity,  16,000  ft/sec) 


The  unsteady  aerodynamic  ablation  due  to  the  thermal  impulse  has  been  treated  by 
Goodman  (7),  who  found  that  total  ablation  will  be  increased  if  the  unsteady  ablation  is 
considered.  Photographs  taken  at  early  stages  of  flight  at  NRL  show  a  luminous  trail  as 
seen  in  Fig.  4.  This  trail  is  caused  by  burning  and  excitation  of  material  which  is  ablated 
from  the  surface.  This  material  collides  with  air  molecules,  and  oxidation  results  (22). 


Fig.  4  -  Luminous  trail  of  an  aluminum 
sphere  during  hypersonic  flight 


There  is  an  indication  that  this  luminous  trail  is  connected  with  the  unsteady  heat-impulse 
phenomenon.  An  estimate  of  time  to  reach  a  quasi-steady  state  of  ablation  during  flight 
in  an  environment  (6)  is  given  by  the  following  order  of  magnitude: 


to.  =  O 

q. 


1.32 

^liq 


y(w), 


.  3  o3 


where 


Suj  2 

^  (1  -  0.252  M/ 


-  0.0175  Mj")  . 


The  subscript  2  refers  to  the  values  immediately  behind  the  normal  shock.  Test  condi¬ 
tions  at  NRL  show  that  during  the  first  3  millisecond?  the  sphere  is  subject  to  unsteady 
ablation  due  to  thermal  Impulse.  This  time  amounts  to  the  order  of  10  percent  of  the 
total  flight  time.  Hence,  the  analysis  presented  may  yield  a  total  ablated  mass  which  is 
leas  than  the  measured  value. 
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SYMBOLS 

Reference  area  for  drj^ 

specific  heat 

Drag  coefficient 

Son  friction  coefficient 

^>ecific  heat  at  constant  pressure 

Diameter  of  the  sphere 

Nondimensional  time 

Latent  heat  of  fusion 

Latent  heat  of  evaporation 

Distance  between  the  original  surface  and  the  gas/liquid  interface  (hi  gas/ 
Uquid  interface) 

Enthalpy 

Energy  absorbed  by  heat  capacity  (c^t) 

Energy  absorbed  by  evaporation 
Radiation  heat  flux 
Thermal  conductivity 
Rear  endpoint  of  the  ablated  layer 
Mach  number 

Molecular  weight  ratio  (gas  to  metal  vapor) 

Mass 

Static  pressure 

Static  pressure  at  the  shoulder  of  the  sphere 
Prandtl  number 
Heat  flux 

Radius  of  the  sphere 


17 


18 


NAVAl  R£5£ARCH  LABORATORr 


Reynolds  number 

r  Radial  coordinate  of  the  body  of  revolution 

Radial  coordinate  of  the  body  for  the  center  of  h  (see  Fig,  1) 
s  Distance  between  the  original  surface  and  the  ablated  surface  (si  solid  surface) 
T  Absolute  temperature 

Tj,  Initial  temperature  of  solid  below  melting  temperature 

Stagnation  temperature 
T,  Recovery  temperature 

Tiub  Sublimation  temperature 

T,  Wall  temperature 

t  Time 

u  Velocity  component  in  x-direction 

V*  Volume 

V  Velocity  component  in  y-direction 

V,  Ablation  velocity  (normal  to  wall) 

X  Coordinate  along  the  surface,  measured  from  the  forward  stagnation  point 
y  Coordinate  normal  to  the  surface 

>3  >3=1-  pVp, 

Stagnation  velocity  gradient  parameter 

/ 

y  Ratio  of  specific  heats 

?  Liquid  layer  thickness  (5  =  s-h) 

M  Dynamic  viscosity 

V  Nondimensional  parameter 

P  Density 

T  Shear  stress 

y  Conduction  parameter 

^  Ratio  of  heat  transfer  with  and  without  mass  injecticn 
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SUBSCRIPT 

s  stagnation  point 

«  Edge  of  the  gas  boundary  layer 

0  Nonablating  body 

r«f  Reference 

•  Free  stream  in  front  of  the  shock 
2  '  Flow  conditions  Immediately  behind  the  normal  shock 
b  Body 

liq  Liquid 

•oi  Solid 

i  Gas-liquid  interface 

SUPERSCRIPT 

Ablation  due  to  melting  and  partial  evaporation 


Appendix  A 

ESTIMATES  OF  THE  VALUES  OF  ^3,  AND  4^ 


In  order  to  estimate  the  values  of  fj,  ^3,  and  which  are  found  in  Eq.  (24), 
the  following  conditions  are  assumed.  An  aluminum  sphere  4  cm  in  diameter  travels 
at  M„  =  15  under  standard  atmospheric  conditions.  At  the  sonic  point  (m^  =  1),  the 
thickness  of  the  liquid  layer  is  i  =  0.001  cm.  For  the  estimate  this  thickness  is  re¬ 
garded  as  constant  over  x .  The  flight  time  is  0,003  sec. 

CALCULATION  OF  AND 

Initially  the  values  of  and  are  computed.  From  Eq.  (24), 

(f)^  2  (f)  “"(r)] 

and 

fa  =  (P,-PJ  • 

The  value  of  x/R  at  the  sonic  point  (ji^  =  1)  from  Eq.  (18)  is 


The  pressure  at  the  stagnation  point  can  be  computed  by  using  normal  shock  relation¬ 
ships.  Then  * 

=  0.561x105  (Ai) 

and 


fjj  =  1. 74x10®.  (A2) 


ESTIMATION  OF  ^  AND  ^ 

3  *4 

From  Eq.  (24) 

=  aR  sinJI-jGjf*) 


(A3) 

(A4) 
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where 


Cjfx)  =  X 


In  order  to  estimate  the  values  of  and  It  is  necessary  to  compute  a  from 
Eq,  (21a); 


The  value  of  f  is  estimated  as  follows: 


^ = m'  -  m- 


Iii5  =  !!i 

r._  c. 


^Uq 

/“liq' 

^ig  ‘ 

where  c  |  is  the  skin  friction  coefficient.  In  laminar  flow 

where  c  Is  a  ccHistant  and  ■I  is  the  length  of  the  gas  boundary  layer  or  liquid  layer. 
Therefore, 

/  t  3/2  1/2  .  1/2 


«g  \  ^**llqN  I  ^ig  V 

^liq)  V“»g/  \^ig/  \^Hq/ 


For  the  estimation  of  f  ,  the  right-hand  terms  of  this  equation  are  computed  as  follows: 

/_  /-»  r_ _ _iTi/2 


"  \T„y  L  288  J 


=  6.21 


From  Eq.  (20) 


Mjiq  =  1.704  centipoise. 


^nce  =  0.018  centipoise, 

(M- 

=  p,.*.,  (1  +  0-2  mJ  =  4/)^. 
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Therefore, 


By  assuming  the  total  length  of  molten  surface  L  to  be  L  =  S  5, 

L  =  (“Hq). =  55 


C 


or,  for  t  =  0.003  sec  and  5  =  O.OOl  cm, 
^  “  1.667  ctn/sec. 

Also 


=  M.a.  =  rf 


.(1  +  0.2  m/) 


■1/2 


Therefore 


=  2. 09  *  10  ■  ®  cm/ sec  . 


/u 

(■C)  ■ 


-9 


For  further  computation  of  average  values,  the  length  of  the  molten  surface  is  assumed 
to  be 


‘'Uq  -  •2  • 


Then 


f  =  1.09*  lO'^K  . 

Because  the  average  value  of  U|j^  is  very  small,  the  liquid  layer  can  be  considered  to 
be  frozen,  so  K  -•  1  and 

f  =  1.09*10-7.  /  (A7) 

Another  term  of  Eq.  (A5)  to  be  evaluated  is 


p,  -  /  /MVT.j 


Since  =  14.9xlO'^cmVsec  at  IS^C, 


(AS) 
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the  assumption  that  the  flow  over  the  sphere  is  laminar,  Newton’s  impact  the  assumption  that  the  flow  over  the  sphere  is  laminar,  Newton’s  impact 

theory  is  applicable  and  the  gas  does  not  dissociate;  the  total  ablated  volume  theory  is  applicable  and  the  gas  does  not  dissociate:  the  total  ablated  volume 

and  shape  of  the  sphere  can  be  computed  as  a  function  of  flight  trajectory  by  and  shape  of  the  sphere  can  be  computed  as  a  function  of  flight  trajectory  by 

using  the  analysis  developed  in  this  report.  using  the  analysis  developed  In  this  report. 


